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Abstract 

We consider the decay of “false kinks,” that is, kinks formed in a scalar field theory with a pair 
of degenerate symmetry-breaking false vacua in H-l dimensions. The true vacuum is symmetric. 
A second scalar field and a peculiar potential are added in order for the kink to be classically 
stable. We find an expression for the decay rate of a false kink. As with any tunneling event, the 
rate is proportional to exp(—S'^;) where Se is the Euclidean action of the bounce describing the 
tunneling event. This factor varies wildly depending on the parameters of the model. Of interest 
is the fact that for certain parameters Se can get arbitrarily small, implying that the kink is only 
barely stable. Thus, while the false vacuum itself may be very long-lived, the presence of kinks can 
give rise to rapid vacuum decay. 
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I. INTRODUCTION 


Kinks are topological solitons in 1+1-dimensional field theories with a real scalar held 0 
and spontaneously-broken discrete symmetry 0 —>■ —<p. The potential has degenerate vacua 
at held values 0 = ±n; the kink interpolates between the two. The same model in higher 
dimensions gives rise to extended objects: linelike defects in 2+1 dimensions, domain walls 
in 3+1 dimensions, and so on. 

We are interested in models which “demote” the vacua at ±n to false vacuum status, 
there being a lower-energy true vacuum at 0 = 0. A specihc example is 

( 1 ) 

where (see Fig. 

Ui(0) = A(02-^n2)(02-n2)2 (2) 

with 0 < (5 < 1. In the true vacuum the symmetry is of course restored. The motivation for 
studying such a model is that the presence of topological defects can have a dramatic ehect 
on the quantum mechanical stability of the false vacuum. Using cosmological language for 
convenience, if the universe is in a false vacuum (without kinks) throughout space, it will 
decay through quantum tunneling [T], yet the decay rate per unit volume can be exceedingly 
small, to the point where the observable universe could be trapped in a false vacuum for 
times exceeding the age of the universe. Such a scenario is invoked in certain models of 
fundamental physics; see for instance [2]. If so, the presence of topological defects can have 
a dramatic effect on the decay rate. (It should be noted that if they exist, topological 
defects will necessarily be formed during a phase transition, so their presence is not merely 
a possibility; it is a certainty |3].) This situation has been examined previously for magnetic 
monopoles |1], vortices in 2+1 dimensions [S], and cosmic strings [B]. 

Since there is a unique true vacuum in the class of model we are considering, there 
are no classically stable nontrivial static solutions. For instance, a kinklike conhguration 
interpolating between the two false vacua would bifurcate into two halves, one interpolating 
between 0 = —v and 0 = 0 and the other between 0 = 0 and 0 = n. The potential energy 
density being lower between the two halves of the kink than in the exterior region, the kinks 
would essentially repel each other and fly off to spatial inhnity at speeds approaching that 
of light, leaving an ever-expanding region of true vacuum. 
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FIG. 1. Example potential Q with symmetric true vacuum and symmetry-breaking false vacua. 
Here A = u = 1, 6 = 0.1. 

Thus we must consider a more complicated model if we wish to study the decay of 
classically stable kinks formed in a symmetry-breaking false vacua. The paper is outlined as 
follows. In Section 2, we introduce a generalized model and argue that it does have classically 
stable kinks. In Section 3, we find numerical solutions for a variety of parameters. These 
solutions, although classically stable, will tunnel to the true vacuum, a process analyzed 
in Section 4. The comparison between kink-mediated vacuum decay and ordinary vacuum 
decay is analyzed in Section 5. We then present conclusions and suggestions for further 
work. 


II. A MODEL WITH CLASSICALLY STABLE KINKS 

As argued above, classically stable kinks do not exist in the simplest model in which the 
vacuum structure is as outlined above (two symmetry-breaking false vacua, symmetric true 
vacuum). One way to obtain stable kinks is to add a second scalar field to the model with an 
unusual potential. Consider the following model, which is to be viewed more as an example 
having the desired vacuum structure and stable kinks rather than as a realistic model for a 
particular physical system. 


I^ = \ 
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where 


^(0,X) = - 5ivi^){(l)^ - vi^f + ^ 2 [(X^ -'1^2)^ - (<52/4)(x - 2t;2)(x + 'i'2)^] 

0 ^ + '-yvi '- 

The model is symmetric under 0 —>■ —0. The potential has seven parameters, two of which 
can be eliminated by rescaling the helds and x to dimensionless variables. Doing so, we 
can rewrite the (dimensionless) Lagrangian as follows, having chosen the constants in such 
a way as to simplify the equations of motion: 

1 ry 

C=-^{d,cj>Y + ^^{d,x?-V{cj>,x) ( 3 ) 

where 


1/(0, x) = (02-hi) (0^-1)^ + 


a 


+ 7 


(X^- 1 )^-^(X- 2 )(X + 1 )^ 


( 4 ) 


There are now hve parameters, a, 0, hi, ^25 7 , which we take to be positive; furthermore, we 
suppose 0 < hi < 1 and 0 < h 2 < 16/3 in order for the potential to have the properties we 
desire. Although it entails a loss of generality, in what follows we will assume a = 0 = 1 for 
simplicity. 

The potential is the sum of two terms. The first term depends on 0 only and is, apart 
from rescaling, the potential ([^ of the original model (see Fig. [^. The second term is a 
product of two factors. The second of these (in square parentheses, written V 2 (x) below; see 
Fig. depends on y only. With 82 in the above-mentioned range, V 2 has two minima: a 
global minimum, of zero energy density, at y = — 1 and a local minimum, of energy density 
h 2 , at X = -|-1. The hrst factor can be viewed as a modulating function which varies the 
“strength” of the second factor depending on the value of 0. In particular, if 0 = 0 (its true 
vacuum), the modulating factor is maximal, so that if y passes from one minimum to the 
other, the cost in potential energy where y — 0 will be large. 

The relevant features of the potential, illustrated in Fig. are as follows. The true 
vacuum of the model is (0, y) = (0, —1); its energy density is 1/(0, —1) = —5i. There are 
false vacua at (0, y) = (±1, —1), of energy density V (±1, —1) = 0. There is also a maximum 
in the vicinity of (0, y) = (0,0) which can be quite pronounced if 7 <C 1. Finally, there are 
possible extrema at y = 1 which are less important but not entirely irrelevant to us, so it is 
worth examining briefly the potential there: 

^2 


1 /( 0 , 1 ) = (02 - 50 ( 02 - 1 )^ + 


+ 7 


( 5 ) 
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FIG. 2. Last factor of the potential Here ^2 = 0.1. 


This is the sum of the potential depicted in Fig. [T] and a Lorentzian function. It is of 
interest to ask if V (0,1) can be negative, since as we will see it gives rise to an instability 
of the soliton (though not the type of instability which is of primary interest to us). If the 
Lorentzian is broad and sufficiently large in amplitude, it will raise the potential near 0 = 0 
so that it is nowhere negative (for y = 1). However if the Lorentzian is narrow, it can raise 
the potential at 0 = 0 to a positive value while leaving negative regions on either side. Thus 
there are two cases. In the hrst case V (0,1) is minimized at 0 = 0 so that (0, x) = (0) 1) is 
a local minimum of the potential. This occurs if 


7 > 


1 + 2di 


( 6 ) 


In this case the minimum of H(0,1) is negative if 


^2 < ^i7- 


(7) 


If (|^ is not satished, V (0,1) is minimized at a pair of nonzero values of 0 straddling 
0 = 0. These points are local minima (false vacua) of the potential while (0, x) = (0,1) is a 
saddle point. Now the minimum of H(0,1) is negative if the following surprisingly unwieldy 
condition is satished: 

9(5i^ 


1 
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FIG. 3. (a) 3-dimensional plot and (b) contour plot of G((/>, x) for (5i,(i2,7) = (0.3, 0.1, 0.7). (These 
parameter values do not produce stable kink solutions; they are used here simply to illustrate 
the form of the potential.) The true vacuum is at (</>, x) = (0,-1); there are false vacua at 
{4>,x) = (=tl, —1)- There are also extrema along the line x = 1) the nature of which depend on 
the values of the parameters, as described in the text. (In the figure, they are local minima - false 
vacua - at {4‘,x) = (0,1) and very near (±1,1).) The potential energy density is maximal, and 
can be very large, near (cp, x) = (0, 0). 

We can understand qualitatively why there could be stable solitons interpolating between 
the two false vacua (0, x) = (=tl) ~1)- Consider a conhguration of the form illustrated in 
Fig. 0 where 0 has half-kinks at one interpolating between 0 = —v and 0 = 0 and the 
other between 0 = 0 and 0 = u, these being “enveloped” by a kink-antikink of x at 
where (To simplify the discussion, we will call all of these objects kinks.) There are 

five regions where the helds are approximately constant, two pairs of which are related by 
symmetry; these regions are denoted (i) (between the two 0 kinks), (ii) (the regions between 
the 0 and x kinks), and (iii) (exterior to the x kinks). In these regions the energy density 
comes entirely from the potential energy and is easily evaluated: 

Cp) = —hi H— S2, Vui) = - — ■ — S2, Viiii) = 0. 

It is easy to see that the configuration depicted in Fig. [^cannot be a solution. Consider 
a family of such conhgurations parameterized by where the positions but not the 
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FIG. 4. Configuration to illustrate intuitively the existence of stable solitons. 

shapes of the transitions vary. For the conhgnration to be a solntion, its energy mnst be 
stationary as a fnnction of For small displacements A/^, A/^, the variation in energy 

is 

AE = 2(V(i) - V(ii))Al^ + 2(V(n) - V(ni))Al^. 

Since Vf^m) < the energy is not stationary in indeed, if left to evolve dynamically, 
the two X kinks wonld move towards the origin to rednce the static energy. Similarly, since 
< V{i) if is energetically advantageons for region (i) to collapse to zero, converting false 
vacnnm to trne. 

This variational argnment is not powerfnl enongh in itself to determine if a conhgnration 
snch as that depicted in Fig. wonld evolve into a stable one, becanse as soon as the kinks 
overlap they interact and the energy is no longer a straightforward fnnction of /^. If the 
conhgnration can be deformed to an nnstable one withont enconntering an energy barrier 
along the way, it will not evolve into a stable one. 

One obvions way to go from Fig. |^to an nnstable conhgnration wonld be to deform a. 
constant, x(a;) = —1 (after which the 0 kinks will hy apart, as described above). This conld 
be done by moving the x kinks towards one another so that they annihilate, or alternatively 
by deforming the valne of x between the kinks from +1 to —1. In either case, the helds pass 
throngh the large potential energy barrier at (0, x) — (0, 0), so there can be a large energy 
barrier preventing either deformation from occnrring. 

Of conrse, this argnment is merely snggestive of the existence of stable solitons. In the 
absence of a convincing analytic argnment for their existence, we mnst resort to nnmerically 
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solving the equations of motion. This is described in the following section. 


III. NUMERICAL SOLUTIONS 


The static equations of motion that follow from ([^ (with a = {3 = 1) are: 


0" - 20(02 - 1)(302 - 1 - 25i) + 


Given a static solution, the energy is 


20 


(02 + 7 )^ 


(X^-1)^-^(X-2)(X + 1)^ 


= 0 (9) 


x" - 


(02 + 7) 


(x2-1)(47-352/4) =0 (10) 


x] = j dx Q (0'2 + x'2) + V (0, x)^ . 


( 11 ) 


We look for kinklike solutions interpolating between the false vacua (0, x) = (±1, —1). We 
expect 0 to be odd and x even under space reflection, so we can solve the equations on the 
half-line a; > 0 with the following boundary conditions: 


0(0) = 0, x'(0) = 0, 

</)(x) —)■ 1, x(^) —t —1 as X 00. 


( 12 ) 

(13) 


The numerical approach used is an adaptation of the relaxation algorithm explained beau¬ 


tifully in [7j. Of course, numerically we do not integrate to inhnity so (13) must be handled 
differently. We integrate to some suitably large x^ax- The most obvious boundary condi¬ 
tion would be to impose (0, x)U=xmax = However, this would produce a solution 

which diverges if extrapolated beyond x^ax- Instead, we linearize the equations ( QpO ) about 
(0) x) = (1) ~1) and insist that the fluctuations tend to zero exponentially as a; —)■ cx). 
This gives linear relations between (0, and between (x, xOa:=a:max which we adopt 

as boundary conditions at Xmax- With these, the algorithm produces a solution which would 
extrapolate to (0, x) = (!)“!) as x —)■ cx). 

Figureillustrates a typical solution, displaying the functions 0(x), x{^) ia (a) and as 
a parametric plot in the 0x plane superimposed on a contour plot of the potential in (b). 
The latter is interesting because the equations of motion ( 9p!0 ) have a mechanical analogy: 
if X is interpreted as a time coordinate and (0, x) as Cartesian spatial coordinates, they are 
the equations of motion of a particle of unit mass moving in a potential —V (0, x)- Thus, for 
instance, the fact that 0 goes beyond 1 at around x = 3 (which at hrst sight might appear 









suspect) is actually perfectly reasonable, since the gradient of V (0, x) along the positive 0- 
axis points towards the origin beyond x = 1. In other words, the particle of the mechanical 
analogy feels a force towards the origin, giving rise to the gently curved trajectory as the 
particle’s position crosses x = 0 (see Fig. [^b)). Also of interest, as we shall see, is the fact 
that in the centre of the soliton y reaches +1 and is essentially constant while 0 sweeps from 
near —1 to +1. 



(a) Kink (5i =0.01, 52=0.01, :)fcO.01) 



FIG. 5. (a) Graph and (b) parametric plot of solution for parameters (5i,52,7) = (0.01,0.01,0.01), 
looking much like that displayed in Fig. |^if the widths of the regions (i) and (ii) depicted in that 
figure were collapsed to zero. 

A search for stable kinks was undertaken over a wide range of values of ^1,52 for six 
values of 7. Where found, as a rule they look much like the one displayed in Fig. The 
kink energy is displayed in Fig. Although of course the details vary from one graph to 
the next, they share several striking features: 

1. The energy is almost independent of ^i, and increases as a function of 82 - 

2. For large (approaching its maximum value of unity) and 82 <C 1, there is no solution; 
this region increases with 7. (The nature of this stable/unstable transition will be 
explained below.) 

3. For 82 somewhere in the neighbourhood of unity (the value depending on 7 but very 
nearly independent of 5i), there is no solution. 
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Stable kink energy {y=0.001) Stable kink energy {y=0.003) 



Stable kink energy {y=0.100) 


Stable kink energy {y=0.316) 




FIG. 6 . Energy as a function of Si, S 2 for six values of 7 . Where no solution was found (for example, 
the near corner of each plot), the energy was set to zero. 

The latter two points are best illustrated in a plot of the stability regions in the 6162 plane, 
for the various values of 7 , shown in Fig. The six lower curves indicate an instability of the 
kink of the type discussed earlier (see the discussion around Eqs. (| 6 ]j^). As one approaches 
the stability boundary (imagine decreasing 62 from above), there is an easily-overlooked 
change of behaviour of 0 at the centre of the soliton. Two possibilities are observed (see 
Fig.ga ,b)), depending on whether or not (§ is satished. If it is, the slope of 0 decreases to 
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Stability boundaries 



FIG. 7. (color online) Stability region (between the curves) in the 6162 plane for each value of 7 
considered. 

zero at the origin, while if it is not, 0 develops nonzero flat regions jnst off-centre (negative 
to the left, positive to the right). In both cases, 0 is settling in towards the minimnm (or the 
pair of minima) of V (0,1). Below the stability line this minimnm is of lower potential energy 
density than the false vacna at spatial inhnity, and the kink becomes nnstable, flying apart 
and leaving a lower-energy false vacnum (or pair of false vacna separated by a kink with 
y = 1 and 0 interpolating between the false vacna) in its wake. In fact, each of the lower 
lines in Fig. [^is somewhat blurry because it is not one line but two that virtually coincide. 
One of these lines is the analytically-calculated stability line mentioned earlier (that is to 
say, (|5 satished and ([^ saturated, or alternatively (|^ not satished and ([^ saturated). The 
other line is the stability line found by the numerical scan of parameter space as described 
above. Obviously the excellent agreement between the two is a pleasing conhrmation that 
our numerical work is behaving as expected. 

The six upper curves indicate an instability of a very different type. As 62 increases, the 
potential 1^2 (y) (see Fig. becomes more and more asymmetric and the barrier between 
y = -|-1 and y = — 1 is reduced. Because of this, y at the centre of the kink no longer reaches 
-1-1 (Fig. [^c)). The instability, then, is towards a conhguration where y no longer encircles 
the less-imposing potential barrier near 0 = y = 0, after which it is energetically preferable 
for y to reduce to —1 for all x, leaving the true vacuum in the centre of the soliton, which 
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(a) Kink near boundary, (6-|,62,7)=(.0014,.00015,.100) 


(b) 


Kink near boundary, (5i,52,y)=(.800,.118,.100) 




(c) Kink near boundary, (5-|,52,r)=(-800,.814,.100) 



X 


FIG. 8. Kinks near the stability boundary, (a) Lower boundary with inequality Q satisfied, (b) 
Lower boundary with inequality Q not satisfied, (c) Upper boundary. 

expands rapidly, converting false vacnnm to trne. 

IV. KINK DECAY VIA TUNNELING 

As we have seen, kinks are fonnd over a wide range of parameters. They are classically 
stable, bnt since they are bnilt ont of the false vacnnm, they (like the false vacnnm itself) 
will decay via qnantnm tnnneling. We will explore the decay of an isolated false kink in 
this section. In principle, we should solve the Euclidean held equations for the instanton 
(bounce) solution; the decay rate is expressible in terms of the action of the instanton, as 
we will see below. 

However, this is a formidable task; even for ordinary vacuum decay the action has only 
been evaluated in a certain limit, where the false vacuum energy density is only slightly 
higher than that of the true vacuum. Then the so-called “thin-wall approximation” is valid, 
and the held theory problem is reduced to a single degree of freedom: the wall radius. 

Here we will make a similar approximation, reducing the helds to a single degree of 
freedom which, at least near the stability boundary, should be a reasonable approximation 
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to the field theory instanton. The latter has the lowest possible action, giving rise to the 
fastest possible decay rate; by making an approximation we will calculate an upper bound 
to the true instanton action which results in a lower bound to the decay rate of the kink. 

Since we are interested in decay to the true vacuum, it is the second type of instability 
discussed in the previous section which is relevant, where the amplitude of the deviation of 
X from its true vacuum value y = — 1 goes to zero. As such, it is of interest to calculate the 
energy of a one-parameter family of static field configurations built out of the kink solution 
, Xk{x)), where the deformation parameter, written h{t), modulates the amplitude of 
Xk while not affecting (pk- 


Xk Xh{x,t) = h{t){xk{x) -Fl) - 1. 


Xh interpolates between the true vacuum at h = 0 and the kink at h = 1. The energy of 
a static deformed configuration is U{h) = E[(l),Xh]- Direct substitution and straightforward 
algebra yields 

U{h) = Xh^ - Yh^ + Zh^ + j dx (14) 

where 




iXk + 1 ) 

<t^l + 



iXk + 1 ) 
^1 + 


Z=^Y- 2X. (15) 


We note that X, Y and Z (which depend both explicitly and implicitly on hi, 62 and 7 ) are 
positive; thus, the potential has a minimum at h = 1 , as indeed it must. 

Fig, I shows U{h), for 7 = 0.01, displayed in two ways, highlighting the effect of hi for 
three values of 62 in (a) and vice versa in (b). The relative unimportance of hi, already noted 
in Fig. is readily seen in both of these. Changing 7 produces only quantitative changes; 
thus for definiteness in the remainder of this section we will (unless stated otherwise) consider 

6i=-f = 0 . 01 . 

We see that in all cases there is an energy barrier separating the kink, at h = 1, from a 
configuration of the same energy with h near zero. This energy barrier must be overcome in 
a tunneling event. Clearly from Fig. [^b), the energy barrier decreases as 82 increases; this 
is as expected since we are approaching the second type of kink instability alluded to earlier. 
As 62 decreases, the energy barrier increases, which might be surprising given that here too 
we are approaching an instability. However, this instability is created by tunneling not to 
the true vacuum but rather to one of the false vacua at x = +1 which arise for sufficiently 
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(a) Kink energy {y=0.010) 


(b) Kink energy (y=0.010) 



FIG. 9. (color online) Kink energy as a function of the deformation parameter h. 


small ^ 2 ; in which case there is no reason to think that it will be particularly easy (in terms 
of the size of the energy barrier) to deform y to the true vacuum. 

The recipe for calculating the amplitude for quantum tunneling, which is related to the 
decay rate of the unstable state, is well-known [1]. For kink tunneling, having reduced 
the two helds to a single degree of freedom h(f), the task is particularly straightforward. 
From the held theory Lagrangian density, we derive an effective Lagrangian for h by the 
substitution 

0 X —>■ Xh(t) = h{t){xk + 1) — 1. 


This yields 


where 


L = -Mh‘^ 
2 


U{h) 


M = dx{xk + 1)" 


(16) 


and U{h) is as above. Thus, the parameter h can be thought of as the coordinate of a 
particle of mass M in a potential U{h). The tunneling amplitude is related to the Euclidean 
action of a particle moving in the potential UE{h) = —{U{h) — U{1)). The relevant solution, 
the bounce, written hs, starts at h = 1 at Euclidean time r = —oo, rolls down towards the 


turnaround point (see Fig. 10) at r = 0, and returns to h = 1 as r —)■ cx). Using standard 
methods, the action of the kink bounce is 

rl 


^B,k = / dhx/U{h) - U{1). 
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Deformed kink energy (5-|=52='y=0.01) 



FIG. 10. Energy as function of deformation parameter h. 


Now, U{h) — U{1) is a quartic polynomial, and as can be seen from Fig. 10 it has fonr 


real roots: the double root h = 1, /i+, and another (negative) root which we will call h_. 
Thus we can write 

U{h) - U{1) = X{h - l)\h - h+){h - h_y, 
in terms of X and Y dehned earlier, 

Y -2X± ^yY{Y - 2X) 


h+ = 


2X 


(Note that it is easy to show that Y > 2X.) So the kink bounce action is 


SB,k = VSMX / dh{l- h)^/{h-h+){h-h_). 
Jh+ 


(17) 


The integral can be evaluated analytically, although the result is not terribly transparent; 


it is: 


^(1 - h+)(l - h_) 


24 


- 2h+h_ + - 4(h+ + h_) + 4 


+ —— ) iK + ^ ^ /' 2 iog f^i + a/I - h-) - log(h+ - h_)) , 


16 


which can also be written in terms oiY = Y/X\ 


— Ja - -Y (SY^ - 12Y + 16) + ^ log 

24 V 2 V J 16 ^ 


4-F + 2a/4 - 


Y(Y -2) 


(18) 


(19) 
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This expression depends on a very complicated way on the parameters. It is displayed in 


Fig. 11 as a function of 82 for three values of (5i and for 7 = 0.01. The action goes to zero as 
82 approaches the upper limit of the stability zone, as expected. Also as expected (see for 
example Fig. [^, the action is virtually independent of 5i. 



0.05 0.1 0.5 1 


FIG. 11. Kink bounce action as a function of <52 for three values of di; in all cases 7 = 0.01. 
The kink decay rate is given by 


Fk = Ake-''«•^ (20) 

where the prefactor Ak is fairly difficult to calculate and is beyond the scope of this paper. 
Nonetheless we will discuss it briefly in the next section. 


V. KINK-MEDIATED VACUUM DECAY? 


In order to determine if and under what circumstances kinks have a significant effect on 


false vacuum decay, we must take into consideration three factors: the kink decay rate ( 20 ), 


the kink density, and the ordinary vacuum decay rate. Let us examine the latter factor. We 
suppose that prior to tunneling the universe is in one of the false vacua (0, x) = (il) “1); for 
definiteness let us choose 0 = +1. The tunneling will reach the true vacuum (0, y) = (0,-1). 
Since for any 0 the potential V (0, y) is minimized at y = —1 (see Fig. [^, the bounce solution 
will be one of constant y. Thus the bounce is determined by a 0-dependent Lagrangian 
obtained by substituting y = — 1 in (|^, giving a rescaled version of ([^ with potential 
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depicted in Fig. [T] 


C* = - SiW - 1 )". 

Interestingly, the Lagrangian (and therefore any physical implications) depends only on 
As discussed in [1], the bounce is the minimum-action solution of the Euclidean equation 
of motion (writing r for the Euclidean time) 

with 0 —)■ 1 at spacetime inhnity. The minimum-energy solution is rotationally symmetric 
[8], so dehning r = + r^, the bounce solution based on the false vacuum 0 b, v(’") is the 

solution of the following, where the prime denotes differentiation with respect to r: 


0" + -0' = ;^Ei(0), 0'(O) = 0, lim 0(r) = 1. 

r acp r—^oo 

If we interpret 0 as the position of a particle and r as the time, this is the equation of motion 


of a particle moving in a potential —10(0) (see Fig. 12) with a time-dependent dissipative 
viscosity term; the particle starts at rest at a position to be determined and reaches 0 = 1 
at time inhnity. The initial point 0+ must have positive potential energy to compensate 



FIG. 12. Potential for a mechanical analogy helpful to understand the bounce (di = 0.1). 


for the dissipation; a continuity argument [T] indicates that there is always a solution. The 
solution appears as a “thin wall” if <C 1, in which case the solution and its action can be 
determined analytically; the latter is 
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The solution and action can be computed numerically for any the analytical and nu¬ 
merical actions are displayed in Fig. We see that the two agree for 5i -C 1, as they 
must. 





FIG. 13. Vacuum bounce action as a function of 5i, calculated analytically using the thin-wall 
approximation (expected to be valid for di <C 1) and numerically. The two agree in the domain of 
validity of the thin-wall approximation. 


The vacuum decay rate per unit length is given by 


Fv/T = A^e 


( 21 ) 


where again the prefactor is difficult to evaluate; we will discuss it briefly below. 


Given the decay rate of a single kink (20) and the decay rate per unit length of the vacuum, 
we can now make a formal statement regarding the importance of kinks for vacuum stability. 
Imagine that the kink density (that is, the average number of kinks per unit length) is p. 
Then in a universe of length L there are pL kinks, and the decay rate (that is, the rate for 
one of the kinks to decay) is pTFk = pLA]^e~^’^''^. We must compare this rate with that for 
a bubble to nucleate in the vacuum (that is, far from any kink). This is, assuming the total 
length occupied by kinks is much smaller than the total length, (Fv/T)L = A^e~^^'''L. We 
see there is a critical density of kinks where these two rates are equal: 


p. — p-(>S’B.v->S’B.k) 

Ap 

This expression may look questionable because the density must have dimension 
whereas the right hand side appears dimensionless. However, the factors Hk, A^ are pro¬ 
portional to ratios of square roots of determinants, and as we will now argue, they have 
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different dimensions. In the case of it is [9] 


A. = 


[SbA 



det(-02-a2 + F/'(0)U) 


- 1/2 


( 22 ) 


Let us focus on the determinant factors. These are determinants of the second variation 
of the Euclidean action about the solution (bounce in the numerator, false vacuum in the 
denominator). The denominator is straightforward since it is essentially the product of 
eigenvalues of a Klein-Gordon operator, whose eigenfunctions are plane waves. The numer¬ 
ator is more complicated since 0 b, v depends on a:, r. Thus the eigenfunctions are no longer 
plane waves. Furthermore, since the location of the bounce can be anywhere in the xt plane, 
the operator has two zero modes: 0x06,v and 0,-0 b,v The full determinant of this operator 


is therefore zero; however the prime in (22) serves to remove these two zero eigenvalues from 
the determinant. Dimensionally, therefore, has dimension (mass)^. 

In the case of Ak, we have formally a similar expression: 

- 1/2 


A\, = 


(SbmA" 

det' (-02 + U"{hA) 

V 2^; 

det{-dl + U''{h)\H=i) 


(23) 


where he is the bounce conhguration of the deformation variable h as described in the 
previous section. One important difference is that there is now only one zero mode to remove 
in the determinant in the numerator (corresponding to time translation of the bounce); this 
is directly related to the change in power of the prefactor [H]. It also tells us that A]^ has 
dimension (mass). 

There is no reason for the ratios of determinants to do anything pathological, so we can 
imagine that they are on the order of a fundamental mass scale of the problem (that of the 
scalar helds; call it m) to the appropriate power. Thus 


Pc 


m- 


S'b,- 


^27rS'B,k 

It is clear that as the parameters - and therefore the actions - vary, the dominant effect on 
pc is due to the exponential factor (indeed, this is why virtually all discussions of tunneling 


ignore the prefactor). Fig. 11 tells us that S'B,k is essentially independent of over a wide 


range, while it depends substantially on 62 - In contrast, S'b,v (Fig. 13) depends only on (5i. 
Examining these hgures in detail, we see that if (5i is small (much less than 1, say) while 
82 is of order 1, the exponential factor will be miniscule and even a very dilute presence of 
kinks will have a dramatic effect on vacuum decay. 


19 















VI. CONCLUSIONS 


We have examined the effect of topological solitons on vacuum decay in a toy model in 
1+1 dimensions with symmetry-breaking false vacua and a symmetry-restoring true vacuum. 
It is far from automatic that such a model would have solitons; for instance, in the simplest 
such model any conhguration interpolating between the false vacua would be unstable, with 
two “half-solitons” repelling each other and leaving true vacuum in their wake. The model 
we study is not particularly realistic but it does show the existence of models with classically 
stable solitons. 

We find that solitons can indeed have an important effect on vacuum stability. This 
is essentially because the energy barrier which makes the soliton classically stable, on the 
one hand, and that which makes the vacuum itself classically stable, on the other, are 
independent. In particular, while maintaining a large barrier between true and false vacua 
(resulting in a very long-lived vacuum), the barrier stabilizing the kink can be made small. 
In this case, the presence of even a very dilute gas of kinks would cause the vacuum to decay 
rapidly via nucleation of a bubble of true vacuum in the vicinity of one of the kinks. 

Of course, any model with solitons in 1+1 dimensions will have domain walls in 3+1 
dimensions. The decay through tunneling of a domain wall in 3+1 dimensions would be 
very different than that of the corresponding kink in 1+1 dimensions: the wall would de¬ 
velop a bulge somewhere and the bulge would then expand, analogous to the corresponding 
phenomenon with stringlike solitons [6]. This situation is currently under investigation. 
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